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Abstract 



We find the renormalisation coefficients of tlie quark field and the flavour 
non-singlet fermion bilinear operators for the domain wall fermion action, in 
the regularisation independent (RI) renormalisation scheme. Our results are 
from a quenched simulation, on a 16^ x 32 lattice, with (3 = 6.0 and an extent 
in the fifth dimension of 16. We also discuss the expected effects of the resid- 
ual chiral symmetry breaking inherent in a domain wall fermion simulation 
with a finite fifth dimension, and study the evidence for both explicit and 
spontaneous chiral symmetry breaking effects in our numerical results. We 
find that the relations between different renormalisation factors predicted by 
chiral symmetry are, to a good approximation, satisfied by our results and 
that systematic effects due to the (low energy) spontaneous chiral symmetry 
breaking and zero-modes can be controlled. Our results are compared against 
the perturbative predictions for both their absolute value and renormalisation 
scale dependence. 



I. INTRODUCTION 

Renormalisation of lattice operators is an essential ingredient needed to deduce physical 
results from numerical simulations. In contrast with the determination of hadronic masses, 
physical matrix elements can be determined only if the normalisation of the appropriate 
lattice operators can be related to that of the corresponding continuum operators, conven- 
tionally specified perturbatively at short distances. In principle, lattice perturbation theory 



*new address: Department of Physics, Ohio State University, Columbus, OH 43210, USA 



1 



may be used to establish this connection. However, lattice perturbation theory converges 
slowly and the expansion parameter, the square of the lattice coupling evaluated at the 
lattice scale, g{a)'^, decreases only as an inverse power of ln(a). This makes systematic im- 
provement of perturbative results essentially impossible. This convergence may be improved 
when, following ideas from continuum perturbation theory [jT|, a renormalised or "boosted" 
coupling rather than the bare coupling is used as an expansion parameter. Even so, con- 
siderable arbitrariness remains, and in general it is extremely difficult to go beyond one loop 
order in such calculations. To overcome these difficulties, Martinelli et.al. 0] have proposed 
a promising non-perturbative renormalisation procedure. This method has been previously 
used to determine renormalisation coefficients for various operators using the Wilson 0-0] 
and staggered actions [§]. The purpose of this work is to study the application of this 
technique to the renormalisation of the quark field and flavour non-singlet fermion bilinear 
operators for the domain wall fermion action. 

Domain wall fermions 0-0] provide an action, that at the expense of introducing a flfth 
dimension, has a low energy theory with excellent chiral properties while at the same time 
preserving exact flavour symmetry. These good chiral properties lead to a suppression of 
the possible dimension flve terms in the long-distance effective Lagrangian implying that 
domain wall fermions deflne a lattice version of QCD which is off-shell improved to (9(a^). 
As we will see, these domain wall off-shell Greens functions show remarkably reduced lattice 
artifacts. A study of operator renormalisation coefficients for this action is useful, both 



because these numbers are needed for use in practical calculations of physical quantities 112 



and because it provides an excellent test of the chiral properties of the domain wall fermion 
action in practical simulations. In fact, we find that domain wall fermions perform quite 
well for non-perturbative renormalisation with negligible contributions from explicit chiral 
symmetry breaking. This finding is in good agreement with recent work on the chiral limit 
of quenched QCD with domain wall fermions ]jT2|,|l3|. 

Careful operator normalisation is especially important for the domain wall fermion 
method. As is reviewed in Section pi A| , the interpolating field conventionally used to 
create and destroy the physical modes is exactly localised in the fifth-dimension on the right 
and left walls. Since the actual physical modes extend somewhat into the fifth dimension, 
the overlap between the interpolating field and the physical modes will be smaller than one. 
This implies a wave function renormalisation factor (Zq) which differs from one even in the 
case of free fields. For the eigenvectors corresponding to the smallest 19 Dirac eigenvalues 
examined in the quenched, /3 = 6.0 calculation of Ref. ]]T^], this overlap typically varies be- 
tween 75 and 85%. Fortunately, the non-perturbative methods employed here ^ precisely 
include these effects. 

We begin in Section |I| with a brief summary of the main issues involved in applying the 
non-perturbative renormalisation method. In Section JT|, we give the domain wall fermion 
action and discuss the Ward-Takahashi identities it obeys. Section |IV| builds on this base to 
constrain the ways in which explicit chiral symmetry breaking terms may enter low energy 
matrix elements calculated using domain wall fermions. In Section [V| we give the details of 
our lattice simulations. Section ^ describes the renormalisation of the quark propagator, 
and in Section |VI1| we introduce the quark bilinears and compute their renormalisation on the 
lattice in the regularisation independent scheme. After removing expected non-perturbative 
pole terms, we look for effects of explicit chiral symmetry breaking and find that they 
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are negligible. In Section [VII]| , we avail ourselves of the axial Ward-Takahashi identity 
again to compute the quark wave function renormalisation from the conserved vector and 
partially conserved axial-vector currents. In Section ^ we calculate the renormalisation of 
the non-conserved, local axial current from a ratio of its hadronic matrix element to the 
hadronic matrix element of the partially conserved axial current and find good agreement 
with the results of Section [VI I| . In Section ^ we convert the renormalisation coefficients 
to renormalisation group invariant quantities by dividing out the renormalisation group 
running. In Section ^ we discuss the calculation of the quark wavefunction renormalisation 
from the propagator. 

After comparing our non-perturbative results with recent perturbative calculations in 
Section [XI I|, we end with our conclusions. The details of the exact conventions and equations 
used for the perturbative running and matching are relegated to appendices. 



II. NON-PERTURBATIVE RENORMALISATION 

In the following the method of non-perturbative renormalisation introduced in Ref. 
will be studied. This method uses a renormalisation scheme that is defined by a set of 
conditions that mandate the renormalised values of the operators of interest between external 
quark states, in a fixed gauge, at large virtualities. As such these conditions may be expressed 
in any regularisation scheme (and so this scheme is known as the regularisation independent 
(RI) scheme). In particular this allows the renormalisation factors to be defined in the 
lattice regularisation, opening the way for renormalisation factors to be directly calculated 
in numerical lattice simulations. 

While calculating renormalisation factors from lattice simulations neatly avoids the need 
to perform analytic calculations using lattice perturbation theory, which are both challenging 
and poorly behaved, doing so introduces several issues that must be considered: 

• Calculating the matrix elements of the operators of interest between external quark 
states requires a fixed gauge to be used. This allows for the appearance of Gribov 
copies, possibly obscuring the required comparison with continuum perturbation the- 



ory where only the trivial copy appears. Earlier studies |14| of the size of Gribov noise 
in the calculation of a gauge invariant normalisation factor as a ratio of two gauge- 
variant amplitudes suggest this may not be an important difficulty for the parameters 
used here. However, in future work, this difficulty can be avoided by taking two steps: 
i) Impose the regularisation invariant normalisation condition in a sufficiently small 
physical volume that non-perturbative effects are suppressed, ii) Begin the Landau 
gauge fixing procedure from a configuration that is in a completely-fixed axial gauge. 
Taking these two precautions will insure that any effects of Gribov copies will be similar 
to other non-perturbative effects and will vanish as the comparison with perturbation 
theory is done at weaker and weaker coupling. 

Numerical simulations are performed with a finite lattice spacing. This provides a 
natural condition, 

(1) 

(Jj 
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over the momenta range for which a direct extraction of continuum quantities is pos- 
sible. 

• As the renormahsation factors are determined in a non-perturbative calculation the 
contributions of propagating mesons, and in particular pseudo-Goldstone bosons, must 
be identified and removed. These effects may be reduced by working at high momenta, 
with a natural condition for the absence of significant deviations being 

Aqcd < bl . (2) 



Taking the last two points together suggests that this technique relies on the existence of a 
"window" of momenta, 

Aqcd < bl < - , (3) 

for which the predictions of continuum perturbation theory should correctly describe the 
form of the lattice data. In practical simulations however, it has been found that the effects 
of deviations due the violations of both these inequalities must be taken into account [P|,p!5|,[7| . 

Fortunately, near either edge of this window, the form of deviations from perturbative 
behaviour may be predicted. In the case of momenta too low, the initial corrections may 
be described by an expansion in terms of momentum-suppressed condensate terms by use 
of the operator product expansion (OPE). In turn, the first corrections to continuum- like 
behaviour may be taken into account in terms of an expansion in the lattice spacing, a. 

Another trivial consequence of the restricted range of momenta available in current lat- 
tice simulations is the need for many phenomenological calculations to be composed of 
continuum perturbation theory calculations at high scales, that are then run down to scales 
accessible on the lattice and combined with the lattice result. As the majority of the existing 
calculations for the continuum perturbative results use renormahsation schemes that may 
only be defined when using dimensional regularisation (such as the MS scheme), perturba- 
tive matching calculations between these schemes and the ones that may be defined in the 
lattice regularisation need to be performed. 



III. DOMAIN WALL FERMIONS 

In this section the domain wall fermion formulation, as used in our simulations, will be 
reviewed. 



A. Action 



The domain wall fermion (DWF) method is a promising new approach to lattice QCD 
introduced in Ref. |T0[, which, at the expense of introducing an extra, discrete, non-gauge 



dimension, provides drastically improved chiral properties at finite lattice spacing while 
preserving exact symmetry under vectorial flavour rotations. This is achieved by using an 
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action in the fifth dimension that is asymmetric between the left-handed and right-handed 
components of the fermion field. Denoting the fifth co-ordinate as s, with 



s G 0,---,L^ - 1, 

the massless action may be written as 

5fermion(m/ = 0) = - J] (-7^.^ (V+ + V 

+ 

with 



(4) 



(5) 



2 V, V+ + M5 



■ x,s 1 



[rff/][c/^ci^]exp(-5, 



gauge 



fermion 



In Eqs. I and ^ '^x,s is the fermionic field, f/^(x) is the gauge field and 

1. 



^ gauge 



p 



-ReTr [Up; 



(6) 



(7) 



with j3 = Q/qq and Qq is the bare lattice coupling. The projectors for the left and right- 
handed spinors are defined as, 



Pl 
Pr 



(8) 



The notation has been used to denote the discrete forward/backward covariant deriva- 
tives: 



(9) 
(10) 



and 9^ represents the corresponding derivative with no gauge term. For the case of the 
derivative in the fifth dimension, df, the domain wall is implemented by giving the derivative 
hard boundaries. For example a one-dimensional acting on a space with four points may 
be written in matrix form as 





\ 


/ 


-1 


1 







( ^ 











-1 


1 



















-1 


1 




V ^x,3 


/ 


V 











-1/ 





(11) 



It should be noted that the action in Eq. |^ is actually the hermitian conjugate of the action 
proposed in Ref. |jTO|. This change was made for practical reasons related to compatibility 
with the existing Wilson operator implementation for the QCDSP machine. 
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In the free theory, for < M5 < 2, the effect of this is to produce a spectrum with one 
hght fermionic mode, with exact chiral symmetry in the Lg ^ 00 hmit, and IQLg — 1 heavy 
modes. The wavefunction of this hght mode has its right-handed component concentrated 
on the wall at s = — 1 and its left-handed component on the wall at s = 0. This light 



fermion mode may be studied by introducing an interpolating operator of the form |16 



Qx = Pi'^xfi + Pr^x,Ls-i (12) 
= '^xfiPR + '^x,l.-iPl . 

The above considerations also naturally lead to the introduction of an explicit mass term to 
the action of the form 

Sierrmonirrif) = Sfcrmionl"!/ = 0) + ^ 771/ gg , (13) 

X 

where is the bare quark mass. In the free case, this leads to a spectrum with one light 
fermion of mass 

Ms (2 - M5) {mj + (1 - Mg)^^) . (14) 

Note that in the oo limit this is proportional to rrif, while for finite Lg there remains 

a residual mass, mj-es, that acts as an additive renormalisation to rrif. 

However, the properties of domain wall fermions in the presence of gauge fields is a 
much more difficult question. In particular while the form of the mass of the light mode is 
expected to be proportional to mj + mres, the dependence of m^es on Lg must be determined. 
Perturbative calculations [[T^-pTI have shown that the existence of the light mode is stable 
to small perturbations and that this mode has all chiral symmetry breaking proportional to 
rrif as Lg ^ 00. These studies also highlight several issues that must be considered when 
undertaking numerical simulations: 

1. The dependence of m^es on Lg may no longer be of the simple exponential form shown 
in Eq. |T|. 

2. M5 undergoes a strong additive renormalisation. This is understandable, as the five 
dimensional problem has no approximate chiral symmetry to protect it. 



Indeed, extensive numerical studies in the quenched approximation ||12|,|13| have shown 
that the Lg dependence of m^-es does not fit a single exponential in the range Lg = 12 —>■ 48 
for lattices with the same lattice spacing (a = 0.520GeV~^) as the results in this paper. For 
Lg = 16, the value used in this work, ?77res was found to be ~ 4MeV in the MS scheme 
at 2GeV [jl2[. The strong additive renormalisation of M5 requires that an input value be 



chosen numerically so that a single light mode forms and that its decay in Lg is as rapid as 
possible. It has been found that for even coarser lattices than used here such a choice can 



be made 
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B. Lattice Ward-Takahashi Identities 



For the purpose of analysing the consequences of the symmetries of the action, it is 
convenient to introduce an extended mass term, M, with flavour structure such that the 
mass term reads 

q^M^qR + qnMqL , (15) 

and so the mass term is invariant under a transformation of the quark fields and the mass 
matrix M of the form 

qi ULqL 
qn Unqn 

M^UrMUI (16) 



Following Ref. |]T6[, on a finite lattice, an exact vector Ward-Takahashi identity may be 
derived by considering transformations of the 5-dimensional fermion field, such that 

5v^.,s = -te:'^a:,sT\ (17) 

where {T"} is the set of hermitian traceless matrices acting on SU{Nf) flavour-space. This 
leads to an exact Ward-Takahashi identity that reads: 

-d^{V;iz)Oix,,---,Xn)) (18) 
+ {q [M, T"] q{z) 0(xi, ■ ■ ■ , x„)) = -^(5'^0(xi, ■ ■ ■ , Xn)) , 

where 

K(^) = IJ2 {^^^+^,s (1 + 7m) Ul^T^<iJ.,s (19) 

^ s 

^x,s (1 7/^) Ux^fjT ^ x+fJ,,S ^ 

For the case of axial transformations the analogous choice is a transformation of the form: 

SA^.,s = -tels'^x,sT\ (20) 

with 

a el 0<s<Lj2 

ex,s I _^a J^J2 <S<Ls. ^ ' 

This leads to a Ward-Takahashi identity of the form 

-d-^{Al{z)0{x,r--,Xn)) (22) 
+ (g{M,r'^}75g(z)0(xi,---,x„)) 

+2{J^^ 0(Xi, ■ ■ ■ , Xn)) = -l{5\0{Xr, ■ ■ ■ , Xn)) , 



where 



A; = ^Y: sign (s - (v&x+,,s (1 + 7.) Ut,,T^vl;,,, (23) 

Jtq = -'^x,l./2-iPlT''^.,l./2 + ^.,l./2PrT''^.,l./2-i ■ (24) 

Therefore, in contrast to the previous case, the axial current is not exactly conserved. This 
is necessary both to provide a mechanism for physical terms due to the U{1)a axial anomaly 
to enter the calculated amplitudes and also to allow for explicit chiral symmetry breaking 
contributions at finite Lg. The situation is analogous to that for Wilson fermions ||23|| , 
where the role of J^^ is played by the chiral variation of the Wilson term, except that the 
contributions from J^^ are expected to tend to zero as — > cxd in the present case |0. The 
form of the contributions from will be further discussed in the next section. 



IV. OPERATOR MIXING AND CHIRAL SYMMETRY 

The major attraction of the domain wall fermion formalism is its ability to decrease the 
size of chiral symmetry breaking by increasing the parameter L^, the distance between the 
two four- dimensional lattice boundaries to which the left and right chiral modes are bound. 
However, it is often impractical or inefficient to choose such a large value of Lg that all chiral 
symmetry breaking effects from mixing between these walls can be neglected. Thus, it is 
important to characterise the effects of this chiral symmetry breaking and in this section we 
will determine how it can effect the low energy physics of lattice QCD. As we will see, this 
can be done as either an expansion in the size of the wall- mixing effects, which for simplicity 
we will denote by 0(e~°'^^) although the exact dependence may be different, and/or as 
an expansion in the lattice spacing a. 

This analysis is easily made by starting with the interpretation of chiral symmetry pro- 
posed by Furman and Shamir |]16[. Here one identifies the full SU{Nf)L ^ SU{Nf)ii chiral 



symmetry of the continuum theory as the independent SU{Nf) rotation of the fermion fields 
defined on the left- and right-hand halves of the five- dimensional lattice: 

^{x,s) ^ Ul'^{x,s) 0<s<Ls/2-1 
^{x, s) ^ Ur^{x, s) LJ2 <s<Lg-l, 

where Ul and Ur are Nf x Nf special unitary matrices belonging to the left and right factors 
of SU{Nf)L SU{Nf){i. From Eq. |1^ it is clear that this transformation will act on the 
four-dimensional quark fields as a standard element of the full chiral symmetry. 



Of course, the transformation in Eq. ^ whose generators are given in Eq. ^ and Eq. [T^, 
cannot be an exact symmetry of the five- dimensional theory as the derivative terms in the 
fifth dimension, taken collectively, couple the left and right hand walls and prevent such 
independent rotations of this single, five- dimensional field. However, in the low energy sector 
of the theory this symmetry can be quite good. The physical, chiral modes which survive 
at low energy are expected to be exponentially bound to the walls with an overlap that is 
suppressed as Lg increases. The higher energy modes which can propagate freely between the 
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walls are all far off-shell with propagators which are necessarily also exponentially suppressed 
at long distances, especially for the large distance Lg. 

In order to characterise the effects of this controlled symmetry breaking that comes 
from communication between the walls, we will generalise somewhat the Dirac domain wall 
fermion operator of Shamir given in Eq. 0. We will introduce a special-unitary, flavour 
matrix VL in the derivative term joining the four- dimensional planes s = Ls/2 — 1 and 
s = Ls/2. Thus, we will modify Eq. ^by adding the term: 

Sn = -Yl {^x,LJ2-iPl - l) *x,L./2 + ^.,lJ2Pr - 1) ^.■,L./2-l} • (26) 

X 

If we include the transformation of the matrix Q, 

n^Unnul, (27) 
this generalised domain wall Dirac operator will now possess exact chiral symmetry. Note, 



a comparison with Eq. ^ shows that transforms "like a mass term" . 

Thus, if we examine this generalised theory that includes the chiral matrix Q, all am- 
plitudes will become functions of fl but will exactly obey the chiral symmetry described by 
Eq. ^ and Eq. Therefore, we need only understand how the matrix Q will enter the 
low energy Green's functions of interest to determine in a precise way the transformation 
properties of the chiral symmetry breaking induced by mixing between the walls. 

To zeroth order in e~"^^, the fermion degrees of freedom will remain bound to the 
walls and propagation from one wall to the other can be neglected. In such circumstances, 
the matrix Q which is introduced at a point mid-way between the walls cannot enter, the 
amplitude will be independent of fl and hence naively invariant under the full SU{Nf)L ® 
SU{Nf)fi chiral symmetry. To the next order, oc e~"'^% we expect phenomena which involve 
a single propagation between s = and s = Lg — l. Thus, the matrix Q should enter linearly 
in such amplitudes. 

An important application of this analysis is to constrain the form of the effective contin- 
uum action which gives amplitudes that agree with those of the domain wall theory through 
a given order in the lattice spacing. To leading order in the lattice spacing, this effective 
Lagrangian has the standard continuum form. The above analysis requires that the mass 
term in this leading order effective Lagrangian must have the form: 

Zmrufi^^ + c {i^n^Pnij + i^nPLij} , (28) 

where c is a constant with the dimensions of mass. Here the field ip represents a conventional 
continuum multiplet of quark fields and all quantities carry their physical dimensions. The 
first piece is the normal chiral symmetry breaking introduced by the input mass rrif. The 
second comes from mixing between the walls and is required by the extended symmetry of 
Eq. and Eq. to be linear in Q. Thus, this induced mass term can occur to first order 
in the mixing between the walls, permitting c oc e~"^=/a. With the conventional choice of 
Shamir, Qa,b = Sa,b, the second term in Eq. ^ reduces to our usual residual mass term with 
am™ ^ 10^^ O . 



In a similar fashion the 0{a) effective Lagrangian will contain a clover term induced 
by mixing between the walls, again to first order in e~°'^^ , since it also has the permitted 
SU{Nf)L ® SU{Nf)R chiral structure: 
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aci {^a^,F^,n^PR^ + ^a^,F^,VtPL^^] . (29) 

where ci cx e'"^' is 0{am^es)- Thus, such a term is suppressed both by the lattice spacing 
and by the smallness of mres- 

If we extend these considerations to 0{a?) terms in the effective Lagrangian, we can 
conclude that a four-fermi operator of the form 

C2a\^^){^^), (30) 

where C2 is a constant, cannot occur to order e""'^". Since this operator will become a chiral 
singlet only when contracted with two powers of the matrix VL or one power of Vt and one 
powers of the mass matrix M of Eq. the coefficient of such an operator must contain a 
double suppression C2 cx e"^""^^ or a further factor of m/. 



V. SIMULATION DETAILS 



In the following discussions much use will be made of the momentum space quark prop- 
agator in Landau gauge. The first step in calculating this quantity is to fix the gauge. We 
implement Landau gauge fixing by iteratively sweeping over all lattice sites, maximising the 
functional 



E [Tr {U,{x) + Ulix] 



(31) 



At each lattice site we determine a gauge transformation matrix, g{x), which increases the 
value of Eq. Rll The maximum is achieved when 



^Tr \b{x)B\:> 



(32) 



is zero, where 



B{x) = A{x) - A\x) - -Tt{A - A^) 



(33) 



A{x)=J2[ulix)+%{x-fi) 



(34) 



and 



Uf^ix) = g{x)U^{x)g\x + /i) 



(35) 



In practice we stop when the quantity in Eq. ^ is smaller than 10~^. 

On this gauge- fixed configuration the quark propagator, S{x,0), from one source, de- 
noted as 0, to all possible sinks is then calculated. A discrete Fourier transform is then 
performed over the sink positions giving. 



S (p, 0) = ^ exp (-ip'^" ■ x) S {x, 0) 



(36) 
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with 

pT = ' (37) 



where /i is one of x,y,z or t and may in principle he in the range ^ — 1. In practice, 
however, only a subset of this range is used. 

Unless otherwise stated all the data that will be presented is from calculations on a 
16^ X 32 X 16 lattice (where the last number refers to the extent of the lattice in the fifth 
dimension). The simulation was performed at /3 = 6.0 with 2000 heatbath sweeps between 
every configuration and with 2000 thermalisation sweeps performed at the outset. In total 
142 configurations were generated. For this lattice size the momentum range was restricted 
to those momenta for which = 0, 1, 2 for fi = x,y,z and rit = 0, 1, 2, 3. Quark propagators 
for 5 bare masses, rrif = 0.01, 0.02, 0.03, 0.04 and 0.05 were calculated all using M5 = 1.8. 

The results will often be quoted against the square of the absolute momentum, where 
this refers to the Euclidean inner product of the momentum defined in Eq. |3^. To be more 
specific 

M'=E*>r' (38) 

where p is dimensionful. 

VI. THE RENORMALISED PROPAGATOR 

Before we treat the more complicated situation of the fermion bilinears it is necessary 
to first consider the renormalisation of the quark propagator. Neglecting, for the moment, 
potential contributions from lattice artifacts the renormalised quark field may be defined as 

1 

Qrenix) = Zg^qo{x) . (39) 

If we similarly introduce a renormalised mass, defined by 

mren = ^'^"^0 , (40) 

where mo represents a generic multiplicatively renormalised bare mass, then the renormalised 
propagator may be written 

ZqSo {p, mo) Uo=mrcn/Z™ • (41) 

Both Zq and are fixed in the RI scheme by requiring that the renormalised propagator 
obey the Euclidean space relations 

lim -ATrf^lf(p)) =1 (42) 



-0 12 
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While Eq. ^ and Eq. ^ seem to give a simple and appealing way to calculate Zq and 
by directly applying them to the lattice propagators, the effect of both lattice artifacts 

and spontaneous chiral symmetry breaking must be considered. 

Lattice actions with explicit chiral symmetry breaking require an additive renormalisa- 

tion of the input mass, which may be taken into account for domain wall fermions by making 

the replacement 

mo^m/ + mres, (44) 

in the equations above. However, the effects of lattice artifacts on the correct definition of 
the renormalised and improved quark field are more complicated. They have been studied in 
Ref. |2i| for Wilson fermions, where it is noted that there are three terms that may mix with 
the definition at 0(a) in the lattice spacing, giving rise to an expression for the improved 
and renormalised quark field of 

gren = (1 + ^qma) {1 + OCg + VTL^^^) + aCNGI^} ^0 , (45) 

where ^ may appear because the gauge is fixed. If such extra terms appear then conditions 
must be found that allow them to be subtracted from the bare quark field before Eq. ^ and 
Eq. ^ may be applied. In the context of simulations using 0(a) improved Wilson action 
at /3 = 6.0 these terms have been found to give significant contributions to the form of 
the propagator [^. In particular was found to be large. However, they all break chiral 
symmetry and so, following the arguments of Section |^ should be suppressed by a factor 
of 0(ami.es) for simulations using domain wall fermions. As such, studying the form of the 
propagator provides an excellent test of the chiral properties of domain wall fermions. 

The effects of spontaneous chiral symmetry breaking on the form of the propagator are 



well known |25,2q]. The most noticeable effect is that the trace of the inverse propagator 



T^Tr (S-^M) = m^en + C7iM + . . . (46) 



picks up an extra contribution, which at lowest order in a power expansion in may be 
described as 

^Tr [S;JXp)) = ^rcn + Ci^^ 

where, at first order in perturbation theory, Ci = Anas/S. Putting Eq. ^ and Eq. ^ 
together, the predicted form for the trace of the lattice quark propagator is 

^Tr(Sr,,\(ap)) =--- + ^C,Z, (47) 

+ ZmZg {anif + arn.cs} 

+ 2 (cNGI^g - C'q) {apf H , 

where terms of C(mcNGi) have been neglected. 

In Figure |l| we plot the left-hand side of Eq. ^ versus {apY for a variety of values for 
mf. As can be seen, for domain wall fermions this quantity approaches a constant value 
for moderately large values of (ap)^. Also, it is encouraging that while at low momenta the 
effects of spontaneous chiral symmetry are visible, there is no evidence to suggest appreciable 
effects from explicit chiral symmetry breaking. In particular, there is no evidence of a large 
additive mass renormalisation. This is visible in Figure 0, which shows the result of a linear 
extrapolation of the data to the point = 0. Figure |^ shows the slope of this extrapolation, 
which from Eq. ^ is ZmZq at large (ap)^. 
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A. Extracting Zq from the Propagator 



The extraction of Zq from the propagator via Eq. ^ is numerically challenging due to the 
need for a discrete derivative to be calculated. A much simpler method 0, is to calculate 

This quantity may then be related to Zq by a perturbative matching calculation performed 
in the continuum IE? . 



On the other hand, the use of Z' q to determine Zq introduces significant O(a^) errors 
through the choice of how the discretised momenta are defined. If we were to replace in 
Eq. ^ the definition of the lattice momentum in Eq. ^ by, for example, 

p^ = -sinp^a, (49) 

(X 

then the resulting Z' q would differ from that given in Eq. ^ by O(a^) because the trace 
includes an explicit factor of Viil^i- 0ii6 can estimate the size of this error by using various 
definitions for the lattice momentum in the analysis. As will be shown later in Sec. 



this uncertainty is roughly 10-20%. In Sections VIII and IX we introduce two methods for 
computing Zq which avoid this uncertainty. 

VII. FLAVOUR NON-SINGLET FERMION BILINEARS 

A. Introduction 

In the following the renormalisation of the flavour non-singlet fermion bilinear operators 
will be considered. To simplify notation explicit quark flavours [u and c?) will be used in the 
following equations. The most general fermion bilinear may be written as MFjC? with 

Ti G {1, 7^, 75, 7/^75, (yp.v\ , (50) 

where i represents whatever indices the gamma matrices have. The renormalised operator 
is defined as 

[^rrf],,, = \nVd\^ . (51) 

The factor Zr is fixed in the RI scheme by defining the unrenormalised, amputated vertex 
function 

nr,o(p, ^) = ^ / d'z d^x, d^X2 6-^^^-^+*^-^ ( [uVd]^ {z) mo(xi) ^^(xs) )amp, (52) 
the corresponding renormalised, amputated vertex function 

nr,ren(P, = ^'^FfiiP, Q) (53) 
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and requiring that 



Ar„ren(P, P)p^=^^^ = Tr(^^r r ) ^^ ( ^ rinri,ren(p, P)^ = 1 • (54) 



Note, the corresponding, unrenormahsed vertex amphtude is defined by 

] 

THE 

so that 



Ar,o = ^ ^^^r.rj Tr ^ rinr,o(p, P)) , (55) 



Ar„ren(p,p) = ^Ar„o(p,p) . (56) 

While Eq. ^ completely defines a procedure for calculating a renormalisation factor for the 
bilinear operator of interest, practically we need to be able to use a renormalisation condition 
that allows us to match to perturbative calculations. In general the value of Ar^fiip^p) has 
contributions from intrinsically non-perturbative effects (such as those due to propagating 
pions) that perturbative calculations do not include. As we are interested in the value of 
the renormalisation factors in the perturbative regime we either apply the renormalisation 
condition at a high enough momenta such that the non-perturbative effects are suppressed, 
or remove such effects from the data and in the following that is what we will do. We 
will reserve the notation Zi, i G {g, S", P, T, A, V^} for the renormalisation factors in the 
perturbative regime. 



B. Za and Zy 

A (partially) conserved current that is normalised in a fashion which is consistent with 
the usual Ward-Takahashi identities will undergo no renormalisation and the corresponding 
Zy will be unity. In particular, for domain wall fermions and the RI renormalisation scheme 
specified by Eqs. ^ and ^ and imposed at high-momentum, /i ^ Aqcd, we expect 

Za = Zv = 1. (57) 

However, on the lattice the (partially) conserved currents are not local and it is frequently 
more convenient to work with their local counterparts. Provided that these are related by a 
chiral transformation one still has 

Za = Zv. (58) 

This does not, however, mean that A^ must equal Ay, and several mechanisms exist for 
splitting them away from each other at low energies. 

Even if there are no significant effects from explicit chiral symmetry breaking, the effects 
of spontaneous chiral symmetry breaking must be taken into account. Consideration of the 
operator product expansion, to lowest order in powers of shows that A^ and Ay may 
get contributions from terms of the form 
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(59) 



and 

Since such terms, by their very nature, stem from chiral symmetry breaking they are not 
constrained to enter and with the same weight. At large momenta these terms are 
suppressed and do not effect the extraction of Za and Zy- 

If the action being used exphcitly breaks chiral symmetry Za and Zy need not be equal, 
but their ratio will still be scale independent. This means that while A^ and Ay need not 
approach one another at high momenta, their ratio should become scale independent for 
large enough momenta. 

While we want to work in the chiral limit for the extraction of the renormalisation 
factors it is also worthwhile to consider what the mass dependence of A^ and Ay should 
be (especially as we wish to extract the chiral limit from data measured in finite mass 
simulations). Requiring that the "generalised" symmetry introduced in Eq. |TB] is satisfied, 
constrains the mass dependence to either be of the same form as Eq. |5^ (a single power of 
mass multiplying something that breaks chiral symmetry - and therefore by the argument 
of the previous paragraph damped with momentum) or proportional to a second or higher 
power of mass. In the latter case any effect on our data should be negligible. 

The above considerations suggest looking at the quantity 

Aa-Av (61) 

This is shown in Figure ^ and, as with the case for the quark propagator, while the effects 
of non-perturbative breaking terms are visible at low momenta, they are damped at higher 
momenta. There is no significant signal for effects from explicit chiral symmetry breaking, 
since A^— Ay is tending to zero. At momenta of interest, there also seems to be no significant 
splitting due to non-perturbative effects with the difference between Aa and Ay being less 
than 1% in the chiral limit at (ap)^ = 0.8 and smaller for momenta above this. This being 
so it is sensible to use the quantity ^ (A^ + Ay) for the extraction of both ZaIZ^ and Zy/Zq 
to increase the statistical accuracy. This is shown in Figure ^ 



C. Zs and Zp 



For a theory with chiral symmetry the RI scheme preserves the well known MS relations 

Zs = Zp (62) 
Z^ = ^. (63) 

If the potential for explicit chiral symmetry breaking is taken into account, these equalities 
cease to be valid, but the quantities Zs/Zp, ZsZm and ZmZp are expected to be scale 
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independent. However, lattice studies using both the Wilson and staggered actions have 
shown that the ratio Ap/A^, which in perturbation theory would be equal to ZgjZp and 
therefore might be expected to be momentum independent up to small corrections, is strongly 
momentum and mass dependent, with the bulk of this dependence arising from Ap. It is 
instructive to consider the source of this discrepancy . We start from the continuum 
axial Ward-Takahashi identity which is derived by taking the axial variation of the quark 
propagator. This reads 



(m„ + rrid) J dz{u{xi) [u'yr,d] {z)d{x2)) = l-,{d{xi)d{x2)) + {u{xi)u{x2))i-, ■ (64) 
Moving to momentum space gives 

(m„ + nid) {u [u'j^d] d) (p, p) = 75 {dd) (p) + {uu) (p)75 , (65) 
which in the m„ — ^ 171^ = m limit gives 

2m{u [M75(i] d){p,p) = {75, 5 (p)} . (66) 
Amputating and tracing with 75 yields 

mAp(p,p) = i-Tr(S-i(p)) . (67) 
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Neglecting all lattice artifact terms except the additive mass renormalisation (which is jus- 



tified by the discussions in Section [V^ , this leads to an approximate expression for Apiatt, 



including the first order contribution of spontaneous chiral symmetry breaking: 

{mf + rnres) Ap,iatt {ap, ap) = h -j—:^CiZg + ZmZg {rrif + rrires) , (68) 

neglecting terms of (9((ap)^). While in the absence of the condensate term this equation 
reduces to Zp = l/Z^, the condensate term, which is clearly visible in Figure |^, gives rise 
to a pole term of the form 

C,Z, (69) 



(m/ + rrires) (ap) 



in Apiatt- Figure || shows the data for Apiatt with this effect clearly visible, with the rise for 
small [ap)"^ becoming much more pronounced as the mass decreases. 

A similar argument may be put forward for A5. In this case it is the Ward-Takahashi 
identity arising from a vector rotation of the fields. 



(m„ — rrid) J dz{u{xi) [ud] {z)d{x2)) = {d{xi)d{x2)) — {u{xi)u{x2)) (70) 
Moving to momentum space, this gives 

(m„ - rrid) {u [ud] d) (p, p) = (dd) (p) - (uu) (p) , (71) 



16 



which in the — hmit tends to 



{u[ud]d){p,p) = -^S{p). (72) 
Finally, amputating and tracing gives 

1 gTr[g(p)-^] 

"^^=12 • ^^^^ 

Note that this relation should be exact for domain wall fermions (with m = rrif) for any 
value of Lg. Using Eq and noting that both the residual mass and the renormalisation 
factors should be independent of mj, gives the approximate expression 

A -7 7^ CiZ. da^qq) .^.^ 
^SM^-ZmZ, + j-^,^^y (74) 

If the mass dependence of (gg), for small masses, is proportional to only positive powers of 
the mass, then the second term in Eq. ^ is almost certainly unimportant as it is suppressed 
in exactly the region of parameter space in which we are working: large momenta and small 
masses. (The effect might be larger than naively expected, as (gg) is quadratically divergent 
in the lattice spacing.) It is necessary, however, to consider the effects of fermionic zero- 
modes on (gg) . Assuming a theory with chiral symmetry, the spectral decomposition of the 
quark propagator leads to an expected form for (gg), on a single configuration C, of 

- {qq)c = Tr[{Ip + m)-^] (75) 

where hq is the number of fermionic zero-modes, V is the four dimensional space-time volume 
and the A„ are such that IpilJn = i^ni^n- The number of such zero modes should grow more 
slowly than the volume, and so the first term in Eq. |7^ will vanish in the infinite volume 
limit. However, for the lattice parameters used for this simulation the effects of zero-modes 



have been found to be noticeable in both (gg) and hadronic spectrum calculations [0 and 



so must be considered for the present case. Comparing Eq. |7^ and Eq. ^ shows that, as 



m ^ for fixed momentum, A5 gets a large contribution from zero-modes of the form 

-Ak)^, (77) 

that must be subtracted before Zs/Zq may be calculated from Eq. |5^. Figure shows our 
data for A^iatt- While the effect of the condensate term is smaller than that for Zp, it is 
noticeable for the lighter masses. 



D. Fitting the Pole Terms 

Considering Eq. |7^ and moving to lattice notation, the method for extracting Zs from 
A5 becomes clear. Working at a fixed momentum, A5 may be fitted to the form 
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As,latt = ^Jff -2 + C2,S + C3,5(am/)2 , (78) 

[amf + U[amres)) 

with Zq/Zs being given by C2,s- While one might naively expect the denominator in the 
above equation to be aruf + amres, as shown in Ref. |jl2| the residual chiral symmetry 



breaking effects that appear in the pole term in qq are not parameterised precisely by amres 
since the singular behavior of the pole enhances what are expected to order variations 
in the quantity arrij-es- One only knows that the residual chiral symmetry breaking effects 
are of 0{amj.f.s)- However, for the pole subtractions in this paper we have used precisely 
arrif + amres- This is justified since the statistical errors on our data are such that the fits 
are insensitive to the exact value of mj-es- 

The situation for the Zp extraction is slightly more complicated. Examining Eq. ^ and 
Eq. |7^, shows that Ap should have a double pole of the form 

Ap,latt = 7 + 7 T + C3.P + C4,p{amff , (79) 

[am f + am,.es) [am j + arrij-es) 

with cs^p being equal to Zg/Zp and the quadratic mass pole due to zero-mode effects in 
{qq). For simplicity, we have again used amf + amj-^s in the first term of the right-hand 
side in Eq. |7^. For practical purposes, however, the need to fit to the quadratic term may 
be avoided by working with m,f > 0.02. Good evidence that the above fitting forms are 
correct is shown in Figure ||. This shows the average, over all the momenta in the range 
0.5 < (ap)^ < 2.0, of the P^r degree of freedom for a correlated fit to the above forms for 
mf dependence with the power of the pole treated as a free parameter. One sees that a single 
pole is favored for the Zp case while a double is identified in the fit for Zs- Further evidence 
is provided by considering the resulting values for Zs/Zq and Zp/Zq. Figure |^ shows a 
comparison between the extracted values of these two quantities. As chiral symmetry would 
predict for Zs/Zq and Zp/Zq, the two quantities coincide at large momenta. This provides 
an excellent test of both the extraction method and the chiral properties of domain wall 
fermions. This can be further seen by comparing Zp and Z^. (as calculated from the trace 



of the inverse propagator), which is shown in Figure |T0|. This product is clearly very close 
to unity. 



E. Z'p 



The tensor density renormalisation, while sometimes neglected in bilinear renormalisa- 
tion coefficient calculations, is a quantity of use to current lattice simulations |^ . An 



extraction of its value will be postponed to Section 0, but for completeness a plot of At is 



shown in Figure 11 



VIII. EXTRACTING Zq FROM THE EXACT WARD IDENTITIES. 

The vector Ward-Takahashi identity of Eq. |1^ is exact at finite lattice spacing. As such, 
the renormalisation coefficient for the conserved vector current defined in Eq. |T9|, is equal 
to unity. Additionally, the considerations of Section |IV| show that through first order in 
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the residual chiral symmetry breaking the extra J"^ in Eq. ^ can be completely absorbed 
into the additive renormalisation of the mass so that the axial Ward-Takahashi identity, for 
low-energy physics, takes on the normal, continuum form. Therefore, the renormalisation 



factor for the axial current defined by Eq. ^ should also be unity to a good approximation. 

The above facts can be used to compute the quark field renormalisation Zq from Eq. |5^, 
as applied the conserved vector and axial currents. For the case of the conserved vector 
current Eq. reads: 



which therefore implies 



Z-iZv^Tr(7^A^)(p)) = 1 (80) 



Z, = -Tr(7'^A^(p)). (81) 



A similar equation also holds for the conserved axial current. 

Because the formulae for the conserved currents contain fermion fields at separate points 
on the lattice as well as summation over all s, the calculation of the matrix elements for these 
operators can be very expensive. We used a random source estimator to compute the part 
of the sum between s = 1 to s = Lg — 2. Also, instead of calculating all four components of 
A^{p) for a given momentum p, we calculate A^{p) for momenta related to p by interchange 
of its 0th component with each of the other three. We then use the equality 

Tr (7'^A^^(p°, p\ p^ p^)) = Tr (7°A° (p°, p\ p^ p^)) 

+Tr(7°A°(p^p^p^p='^ 
+Tr(7°A°(p^p^p^p3^ 

+Tr(7°AO(p^p^p^pO)) , (82) 

to obtain the needed result. As the time to Fourier transform a matrix element is negligible 
in comparison with the calculation of the matrix element itself, this allows us to obtain the 
result with only a quarter of running time of direct calculation of Ay. 

Because the volume 16^ x 32 used in the simulations is not symmetric, strictly speaking 
the above equation is not exact, as the third component of momentum is not related by 
symmetry to the first and second ones. However, this difference is suppressed by two powers 
of the lattice spacing and, in practice, the results obtained for the last three terms in Eq. ^ 
all agree within statistical errors for the momenta used. We used momenta with the first 
three integer components no larger than 2 and the last component equal to or 2. The two 
exceptions were the momenta with integer components (2, 2, 2, 0) and (2, 2, 2, 2) that were 
excluded since they would require usage of momenta (0, 2, 2, 4) and (1, 2, 2, 4). 

Figure |12| shows the difference between Zg calculated using the axial and vector Ward- 



Takahashi identities, while Figure 13 shows the average. As can be seen from Figure |T^, 
while for low momenta the two methods give different results, this difference is damped 
at large momenta as would be expected if this effect stems from spontaneous (rather than 
explicit) chiral symmetry breaking. Again, this provides strong evidence that the effects of 
explicit chiral symmetry breaking are negligible in these calculations. 
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IX. Za from hadronic matrix elements 



As mentioned in Section |V111|, to a good approximation the axial current defined by 



Eq. |23| is conserved, and therefore has a renormahsation coefficient equal to the identity. 
This provides a simple way to calculate the renormahsation coefficient of the local axial 
current operator, Za-, directly from hadronic matrix element calculations. One method, 



that has been used in Ref. ||T2[, is to note that the matrix element of any operator with 
the renormalised axial current is a well defined quantity, which will be independent of the 
interpolating operator for the axial current at distances above the scale of the lattice spacing. 
Therefore, 

Tl^qit^)) = Za{Ao{U) qi^q{t2)) . (83) 



for |ti — t2|ct ^ 1- A full discussion of this method and the results are given in Ref. |]12|, but 
it is useful to summarise the results here. Table | collects together values for Za for several 
Ls values on a (3 = 6.0, 16^ x 32 lattice with M5 = 1.8. The quoted errors are statistical 
only. 



X. RENORMALISATION GROUP BEHAVIOUR 

The previous sections have provided an extraction of the renormahsation coefficients of 
interest taking into account the possible effects stemming from chiral symmetry breaking 
(either explicit or spontaneous). In general, perturbation theory predicts that these coeffi- 
cients may be logarithmically dependent on the momentum scale. Lattice artifacts may also 
cause the result to depend on the definition of momentum on the lattice. When small these 
lattice artifacts will be manifest in the data as added terms proportional to {apf' ■ 

The simplest approach to using the renormahsation coefficients calculated in this paper 
is to take the (mass-pole subtracted) value for Aj {ap) as {Zg/Zi) (ap). However, we need to 
address two significant uncertainties: 

First, this choice assumes, without verification, that the 0{{ap)^) lattice artifacts are 
small. We should attempt to understand the momentum dependence of the amplitude 
Aj [ap) to determine the size of the (9((ap)^) contamination and, if possible, remove it. 

Second, most operators of interest in lattice calculations are ultimately defined using 
continuum perturbation theory. Making this connection requires the use of perturbation 
theory to connect the RI scheme and momentum scale used here with the renormahsation 
procedure and momentum scale used in the original continuum definition, typically the MS 
scheme. This final matching step between the normalised lattice and continuum operators 
is done at a specific momentum scale for the renormalised lattice operator. In general, 
both the normalisation of lattice operators and the matching coefficients will depend on this 
momentum scale. It is not known a priori how many loops in perturbation theory must be 
calculated to correctly describe the momentum range probed in current lattice calculations, 
or even if perturbation theory can describe the region we are studying. Comparing the 
momentum behaviour predicted from perturbation theory to that of the data therefore 
provides an important consistency check for the general framework of the method. 
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Our approach to comparing the known perturbative running of the quantities of interest 
to our numerical data will be to divide the data by the predicted renormalisation group 
running, with the overall normalisation set by requiring that at the point {apY = 1 this 
divisor is one. If the perturbative result correctly describes the data, and the effect of (ap)^ 
terms may be neglected, the result will be completely scale independent. 

There are three components that are needed to calculate these quantities: 

1. The anomalous dimensions for the operators from perturbation theory. 

2. The ratio of Z^^ /Z^^ must be known to perform the matching. Since the renormal- 
isation condition that determines Z^^ is well defined both on the lattice and when 
in continuum dimensional regularisation this ratio may be calculated perturbatively 
using the latter regularisation. In general it can be expanded as 

£1 = 1 + ^4^)^^ + + . . . . (84) 

For a consistent treatment this ratio need only be known to one less power of as than 
the running is known. 

3. A lattice value for as- The value of as affects the scale dependence of both the matching 
and running for this calculation. For this work the value of as was calculated at three 



loops using a lattice value of Aqcd taken from Ref. as 



Aqcd = 238 ± 19MeV . (85) 



To do this consistently with the way the lattice treatment in Ref. was performed, 
their value of ro = 0.5 fm was taken and converted into a lattice spacing using the 
results of Ref. [^. For the dimensionful scales that we will quote, we set the physical 



scale through the rho mass computed with domain wall fermions [Q, which for /3 = 6.0 
gives 

a = 0.520(ll)GeV~^ (86) 



Both Zyi and Zy should be scale independent, but this is not the case for Zg. Figure 14 
shows both 1/Aa and the scale invariant (SI) quantity calculated as described above: 

A'J {(apf) = Aa {(apf) /Ca {{apf) . (87) 

The quantity Ca is determined through three loops using the anomalous dimension coef- 
ficients calculated in Refs. []32|,^,^ as described in Appendix |B[ It is normalised so that 
C^(l) = 1. As can be seen, in this case the renormalisation group running actually goes in 
the opposite direction from the data. The scale dependence of this data, either predicted or 
actual, is, however, very small and a plausible explanation for this is an {apY error. Indeed, 
when a linear fit of the SI data versus {apY is performed, for 0.8 < {apY < 2.0, the gradient 
is ~ —0.02. A more compelling test of the renormalisation group behaviour is provided by 
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studying the data for Zs/Zq. In this case the predicted scahng behaviour over the range 
of momenta studied is much larger and, as Figure |15] shows, the agreement between the 
predicted behaviour and the data is impressive (with a gradient, in this case, of ~ —0.003). 
The values for Zs/Zq versus momentum used here are taken after the mass-pole has been 
subtracted and, again, the three loop results for the running are taken from Refs. p2|,^,|27 . 



Unfortunately, a matching calculation for Zt could not be found in the literature, so the 
data could only be compared to the one loop running (which is taken from Ref. [Q). The 
SI quantity so calculated is shown in Figure |16| and has a gradient of ~ —0.02. 

Taking the interpretation that the remaining scale dependence is due to 0{{apY) effects, 
the correct way to extract the renormalisation coefficients is to first construct the SI quantity 
as described above, and then fit any remaining scale dependence [Q to the form 

y = ci + C2{apf , (88) 

for a range of momenta that is chosen to be "above" the region for which condensate effects 
are deemed to be important. Table || shows the fitted values for the RI and MS scheme 
renormalisation coefficients using a fitting rage of 0.8 < {apY < 2.0. Now that the renormal- 
isation group running has been taken into account, it is possible to make a comparison of the 
various methods of calculating Zq and thus give final results for the renormalisation factors. 
Table ^ already gives Zq as calculated from the conserved currents (Figure |1^ shows the 



momentum dependence of both the SI and bare form). Another simple way to derive this 
quantity is by taking Z^jZq from Table |I| and combining it with the value of Z^ obtained 
from hadronic matrix elements. This gives Zq = 0.805(3) (15). This, approximately 5%, 
difference may be taken as an indication of the size of the systematic errors. 

XI. Zq from the propagator results 

Here we show results for the wave function renormalisation computed through Eq. ^ 
and demonstrate that this methods contains a comparatively large systematic uncertainty 
due to the ambiguity in defining discrete momentum. We use the perturbative matching 
between Z'q and Zq, as given in Ref. Then the SI Zq can be constructed as described 

above. 



Figure |18| shows the SI Zq using Z'q defined in Eq. and Figure [T^ shows the SI Zq 
where the replacement 

ap^ sin(ap^) = ap^ (89) 

is made in Eq. Note that the former is plotted vs. {apY and the latter vs. {apy. We use 
the data at rrif = 0.02 since no mass dependence can be observed. 

As in the previous section, we extrapolate to (ap)"^ = 0. We find for the data in Figure 



Zq = 0.715 ± 0.007 ± 0.040, where the first error is statistical and the second comes from 
different choices for the range of momenta over which to fit. The data in Figure ^ give 



^Their convention is that a given Z-factor in Ref. |27| is the reciprocal of ours. 
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Zq = 0.733 ± 0.007 ± 0.050. We can further probe these discretisation uncertainties by 
extrapolating to zero (ap)"^ or (ap)^. This results in Zg = 0.732 ± 0.006 ± 0.020 and 
Zq = 0.772 ± 0.006 ± 0.020, respectively. 

The spread in values of Zg obtained depending on momentum ranges and on the definition 
of discrete momenta mean that extracting a SI Zg in the same manner as the Z's for bilinear 
operators is less precise. However, the results are in rough agreement with the more precise 
methods described above. 



XII. COMPARISON WITH PERTURBATION THEORY 

All the renormalisation factors considered above have also been calculated in lattice 
perturbation theory, at the one loop level, for the domain wall fermion action in the Lg ^ oo 
limit |]20| , p!9| , ^ . As we see little evidence of explicit chiral symmetry breaking effects in our 
study, the fact that the perturbative calculations have been performed in the Lg —>■ oo limit 
will probably not affect this comparison. However, a more serious issue is which M5 value 
to use in the perturbative formulae. 

The reason for this is easy to understand. Away from the walls, the massless domain 
wall fermion Lagrangian, Eq. |^, may be viewed as a simple extension of the standard Wilson 
fermion action to five dimensions with a negative mass term, M5. Like the four dimensional 
Wilson mass term, M5 undergoes a strong additive renormalisation, the size of which per- 
turbation theory is not good at predicting. While a more accurate prediction may be made 
using tadpole improved perturbation theory |@], a good deal of ambiguity remains in the 
perturbative prediction of any quantity that is strongly dependent on M5. 

Further consideration of the similarity of the domain wall and standard Wilson actions, 
leads to a non-perturbative estimate of the magnitude of this additive renormalisation. As 
argued in Ref. in the — 00 limit, the effect of this additive renormalisation may be 



taken into account by using 



Mi"'' = Mf"" - ( 4 - — ) , (90) 



2 k 



in the perturbative equations, where Kc is the four-dimensional critical Wilson hopping 
parameter and M|™ refers to the value used in the non-perturbative simulation (in our case 
M|™=1.8). For f] = 6.0, this ansatz leads to shift of magnitude ~ 0.8, which has been found 
to describe the dependence of the pion mass squared as a function of M5 to a good degree 
of accuracy in a numerical simulation with Lg = 14 pl[]. 



Eq. 4.10 in Ref. [1£] gives the complete one loop bilinear renormalisation constants in 



the MS scheme: 



Z'^°'^' = {{l-wl)Z^)-'Zr. (91) 



Here Z^, Z2, and Z^ are to be computed from Eqs. 3.30, 3.42, and 4.11 and Tables 2 and 
3 in Ref. |[T^, while Wq = 1 — M5. In the mean-field improved case, the above relations 
hold with the replacements Wq Wq'^ = Wq + 4(1 — u), Zyj Z^^^ , Z^ uZ^^^ , and 
Zy uZ^^^ whose values can also be computed from Tables 2 and 3 in Ref. [|T^. The 
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factor u in these formulae is the mean hnk variable in Feynman gauge. As it is not possible 
to use the value of the mean link in Feyman gauge we have instead used the fourth root of 
the plaquette and the perturbative results of Ref. to convert the results of Ref. [1^ . 



In Figure |20| and Figure ^ we plot Z^^^^^ and Z^"**^', respectively, as functions of the 
variable M5 in naive perturbation theory, in naive perturbation theory with the variable M5 
shifted according to Eq. ^and in the mean-field improved case. To compute a^, we used the 
same input values for Aqcd and a as in the perturbative running calculations in Section 
We obtain {{apY = 1) = 0.20. These figures show appreciable M5 dependence. Our non- 
perturbative result is shown as a point corresponding to the single value of M5 = 1.8 that 
we have studied. 

The naive perturbation theory curve has a significant dependence on the precise value 
of as- In the mean-field improved case this problem is not as serious as the coefficient of 
as is a factor of 2-3 times smaller. Examining Figures ^ and one recognises that naive 
perturbation theory does a poor job of determining or Zs giving values nearly 2 times 
too small for M5 = 1.8. Introducing the shift of Eq. Q improves the situation noticeably 
giving values 15% too small and to within a few percent, although the perturbative result is 
rapidly varying with M5 in this case. The mean-field results differ from the non-perturbative 
result by around 5% in both cases. 



XIII. CONCLUSIONS 

In this paper we have described a first study of non-perturbative renormalisation of the 
quark field and fiavour non-singlet fermion bilinear operators in the context of domain wall 
fermions. We presented a theoretical argument constraining the form that explicit chiral 
symmetry breaking effects may take, and found that numerically these are insignificant, as 



might be expected from the measured size of the additive mass renormalisation, m^es, |lT2| , p!3 
However, systematic effects due to spontaneous chiral symmetry breaking and zero-modes 
are significant, but accurately follow the expected form and can be effectively subtracted 
away. 

Renormalisation group invariant quantities were obtained in Section ^ by dividing the 
regular isat ion independent scheme coefficients by the three loop renormalisation group run- 
ning (where available). The residual scale dependence of these quantities is small and was 
treated as an O(a^) error. Three different quantities were used to determine the quark renor- 
malisation factor: the off-shell vertex functions of the conserved vector and axial currents; 
the trace of the product of ^^^7^ and the off-shell quark propagator; and the combination of 
Za as determined from hadronic matrix elements with the value of Za/Z^ obtained in this 
study from the off-shell, axial vector vertex function. The technique of obtaining this Zq 
directly from the propagator suffers from large discretisation errors, but is roughly consistent 
with the other two methods which gave results differing by ~ 5%. 

In the final section we compared our results against the predictions of both standard and 
mean-field improved one loop perturbation theory. 
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APPENDIX A: THE RUNNING OF as 

In the following the definitions 



a 



r 

An 



and 



2N, 



c 



(Al) 



(A2) 



will be used. The renormalised coupling may be defined in terms of the bare coupling by, 

ab = Z^gasi^^' . (A3) 
As ctfe is completely independent of /z. 



-eas + P 



with, 



P = -as 



-/i 



,dZg 

djj? 



The results for the beta function are most easily given in terms of the j3i variable: 



An 



-/3o 



'as' 


2 


'as' 


An. 


-A 


An_ 



(A4) 



(A5) 



(A6) 



The values used in the current work are summarised in Table 0. They are taken from 



Ref. with the number of flavours set to zero (as we are working in the quenched case) 
and the number of colours set to three. 

Once the jSi^s are known the running equation. 



(A7) 



may be solved: 
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One loop solution: 

«£_ _ 1 

47r " (3o In (/iVA^ci.) 

Two loop solution p5[ : 

a, _ 1 /3i In In (/iVA^cp ) 

Three loop solution 

as _ 1 Alnln(/iVA^cD 



(A8) 



(A9) 



47r " /5o In (/xVA^cp) f^'o In' (/iVAlcp) ^"^^^^ 

+ /?o^ln3(;.VA^,,) (^'/^^^^) - (^'/^^^^) 

Wo-/??} 

APPENDIX B: THE RUNNING OF THE Z-FACTORS 

As mention previously, the renormalised operators we are working with are defined as 

ZoObare = Or en ■ (Bl) 

Requiring that the bare operator is independent of the renormalisation scale gives the RG 
equation, 



2 



^^Oren- (B3) 



Writing the solution to this equation as 



and using the notation 



7o = i:of(£j (B6) 
7o. = ii- (B6) 



2A 







A = f , (B7) 

gives rise to solutions to the running equation of the form (where we have suppressed the 
subscripts identifying the particular operator O): 
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One loop solution p5[]: 



C(^2)=«,(/ifo. 



(B8) 



Two loop solution: 



C(/.^) = a.(/i)-o|l + ^(7i-^i7o 



(B9) 



Three loop solution: 



(BIO) 



+ 



1 / asifx) ' 

2 I An 



2 , 



_ — _ \ 2 _ — 2_ — _ — _ ' 

7i-/5i7o + 72 + /5i7o - /?i7i - /527o 



Tables to |VI1] show the anomalous dimensions used in this work. These values were taken 
from Refs. |3^^,Q with the number of flavours set to zero and the number of colours to 
three. 



APPENDIX C: MATCHING COEFFICIENTS 



The numerical values of the matching coefficients, Zq^^ and Zq^' in Eq. used for Z, 



A2) 



and Zs are collected together in Table |VII1| and Table |IX. 
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TABLES 



Ls Za no. configs. 

T2 0.7560(3) 56 

16 0.7555(3) 56 

24 0.7542(3) 56 

32 0.7535(3) 72 

48 0.7533(3) 64 



TABLE I. Za computed from the ratios of hadronic matrix elements. 



Z - factor RI/SI 175 Ml 2G(:V 

ZA/Zq 0.934 (2)(10) 0.938 (2)(12) 

Zs/Zq 0.683 (7)(30) 0.779 (8)(35) 

Zr/Zq 1.034 (3) (100) 1.035 (3) (100) 

^Hadronic 0.808 (3)(15) 0.805 (3)(17) 

Zf^'"^ 0.753 (16) (30) 0.750 (15) (30) 

TABLE IL Final Z-factor results. Zq is calculated two ways: from using ZA/Zq from this table 

combined with Za from hardonic matrix elements, denoted ^Hadromc^ from the conserved 
currents using off-shell quark states, denoted ZY'^'^'^. 
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0.7339(42) 
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TABLE III. 


This table collects together the raw data used for the 


"Bare" data plotted in 


Figures 14,15 and 16 for 0.5 < (ap)^ < 2.0. 






Pi Quenched Value 


/3o 






11 


Pi 






102 


P-i 






1128.5 




TABLE IV. Pi's 


for the Quenched theory 




Elements of 7^ 






Quenched Value 

















44.6667 








1056.65 



TABLE V. Quenched Zq Anomalous Dimensions 
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Elements of 75 




Quenched Value 






-8 

-134.667 
-2498 


TABLE VI. Quenched Zs Anomalous Dimensions 


Elements of jt 




Quenched Value 






2.66667 


TABLE VIL Quenched Zt Anomalous Dimension 


7(1) 
7(2) 






-14.4975 


TABLE VIIL Zg RI to MS matching coefficients 


7(1) 

7(2) 
^0 




5.33333 
188.651 



TABLE IX. Zs RI to MS matching coefficients 
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FIG. 1. A plot of ^Tr 

i^^iatt) versus (ap)^ showing that for moderate values of {ap)"^ the effects 
of explicit chiral symmetry breaking are small. 
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FIG. 2. The value of ^Tr {S ^) extrapolated to m/ = vs (ap)^. For moderate (ap)^ the 
extrapolated value is zero within errors, showing that the residual mass is small. 
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(ap) 

FIG. 3. ZjnZq calculated from the slope of ^Tr {S^^) versus m/ plotted as a function of (ap)^. 
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FIG. 4. A plot of — Ay versus {apf', showing that there is no significant difference between 
Za and Zy, even for moderate values of {apf'. 
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FIG. 5. A graph of i {A^ + Ay} versus (ap)^, which up to lattice artifacts, gives Z^jZq and 



Zy/Zq. 
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(ap)' 



FIG. 6. Ap versus (ap)^ for several values of mj, showing that the l/p^ pole is more pronounced 
for small nif. 
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FIG. 7. A plot of A5^iatt versus (ap)"^ for several masses. The mass pole can be clearly seen for 
small momenta and is attributable to zero-mode effects. 
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FIG. 8. This figure displays the average P^r degree of freedom for the fits used to deter- 
mine the power of the mass poles in and Ap, clearly showing their double and single poles, 
respectively. 
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FIG. 9. A comparison of Zs/Zg and ZpjZq as extracted from A5 and Ap after pole subtraction. 
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FIG. 10. The product ZpZ^ calculated by combining Zp/Zq from pole subtraction with the 
trace of the inverse propagator. This product is clearly unity within errors. 
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FIG. 11. This is a graph of At versus (ap)"^, from which Zt wiU be extracted in Section 
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FIG. 12. The difference between Zq as extracted from the conserved axial and vector currents. 
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FIG. 13. The average of Zg from the conserved axial and vector currents. 
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(ap)' 



FIG. 14. A plot showing the raw data for 1/Ka (labeled as "Bare") and the value of 1/A^ 
divided by its predicted three loop perturbative running (labeled as "SI"), such that they coincide 
at (ap)^ = 1, versus momentum. The slope of the latter versus (ap)^ may be interpreted as an 
0{a^) effect and is « -0.02. 



46 




(ap)' 



FIG. 15. A plot showing the raw data for 1/As (labeled as "Bare") and the value of I/A5 
divided by its predicted three loop perturbative running (labeled as "SI"), such that they coincide 
at (ap)'^ = 1, versus momentum. The slope of the latter versus (ap)^ may be interpreted as an 
0(a2) effect and is « -0.003. 
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FIG. 16. A plot showing the raw data for l/Ar (labeled as "Bare") and the value of 1/At 
divided by its predicted one loop perturbative running (labeled as "SI"), such that they coincide 
at (ap)'^ = 1, versus momentum. The slope of the latter versus (ap)^ may be interpreted as an 
0(a2) effect and is « -0.02. 
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FIG. 17. Bare and scale invariant (SI) versions of Zq determined from the conserved axial and 
vector currents. 
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FIG. 18. A scale invariant (SI) version of determined from Z' 
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FIG. 19. The scale invariant (SI) version of Zq determined from Z' ^ but using ap^ instead of 
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FIG. 20. The renormalisation factor Z^^^^ of Eq. |9^ in the MS scheme at 2 GeV computed in 
naive perturbation theory, naive perturbation theory shifted by (4 — 1/2kc) as in Eq. g and mean 
field improved perturbation theory. 
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FIG. 21. Same as Figure |20| but for 
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